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Abstract 

Coherent state transfer is an important requirement in the construction of 
quantum computer hardware. The state transfer can be realized by linear 
next-neighbour-coupled hnite chains. Starting from the commensurability 
of chain eigenvalues as the general condition of periodic dynamics, we hud 
chains that support full periodic state revivals. For short chains, exact so¬ 
lutions are found analytically by solving the inverse eigenvalue problem to 
obtain the coupling coefficients between chain elements. We apply the solu¬ 
tions to design optical waveguide arrays and perform numerical simulations 
of light propagation thorough realistic waveguide structures. Applications of 
the presented method to the realization of a parallel bus for quantum states, 
perfect transfer and edge oscillations are proposed and considered in detail. 
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1. Introduction 

Coherent transfer of quantum states is one of basic requirements for the 
construction of quantum computer [T]. Coherence preserves relative phases 
between quantum states allowing for their interference. The transfer of quan¬ 
tum states can be realised either via moving particles that act as qubit 
carriers, e.g. ions [2], or via stationary solid-state media with qubits that 
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propagate through them, e.g. integrated optical circuits [3]. Scalability and 
portability of future computers have steered investigations towards the lat¬ 
ter, in which qubits are transferred by chains of linearly conpled wavegnide 
arrays (WGAs)[l] or qnantnm dots [5]. Propagation throngh linearly con- 
pled chains is, in general, qnasi-periodic withont fnll revivals of the initial 
state and hence withont a possibility of faithfnl transfer of a qnbit across 
the chain. In isolated noiseless stationary chains, snch dephasing is a conse- 
qnence of incommensnrability of the eigenfreqnencies of the chain. Althongh 
it is a reversible nnitary process and not irreversible decoherence, it prevents 
state reconstrnction and perfect transfer. 

It has been snggested [B] and demonstrated in cold atoms [7j that the 
dephasing can be mitigated by a dynamic control of the propagation along the 
chain by interventions to the system Hamiltonian. However, implementation 
of this method in miniatnre solid state systems that hold a great promise as 
the fntnre qnantnm compnter hardware (e.g. integrated photonic circnits, 
qnantnm dots) reqnires hne interventions along the chain and represents a 
major technological challenge. Alternatively, conplings between the elements 
of a translationally invariant chain can be tailored while leaving the system 
Hamiltonian constant dnring the state evolntion. For example, coherent 
dynamics of qnantnm states throngh a Heisenberg chain of linearly conpled 
atomic Zeeman states [H] and along an invariant optical WGA [1] have been 
engineered by controlling the strength and ratios of their conpling coefficients, 
respectively. Unlike in atomic and other spin systems, in which only the next- 
neighbonr conpling is allowed by selection rnles, in chains withont selection 
rules (WGAs), higher order couplings are present. Since the higher-order 
conpling does not canse decoherence, the periodic dynamics are possible bnt 
reqnire intricate engineering of conpling coefficients [B]. 

Mathematically, conditions of periodic transport and perfect revivals are 
satished when the conpling matrix eigenvalnes are commensnrate mm- 
Hereafter we refer to these chains as commensnrate chains. The chains with 
less than fonr elements are always commensnrate. The eigenvalnes of longer 
hnite chains are commensnrate only for certain ratios of their conpling co¬ 
efficients. The key challenge in chain engineering is hnding these ratios and 
that reqnires a solntion to the non-trivial inverse eigenvalne problem. In¬ 
verse eigenvalne problems are analytically solvable for a small nnmber of 
cases and are, in general, of polynomial complexity |T2]. Inverse solntions 
are not nnique and a concrete choice of solntion depends on the targeted 
transfer dynamics. 
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Realizations of coherent transport in linear chains mostly aim at the per¬ 
fect transfer of the state amplitude (in literature referred to as perfect trans¬ 
fer) without restrictions on phase. As a result of extensive theoretical studies 
of the perfect transfer, a general solution [13] and several particular solutions 
have been proposed mum. Curiously, reported practical realizations mostly 
rely on solutions found in nature in the form of Bloch oscillations [12] , atomic 
spin chains [8], or their emulation by optical lattices [16]. However, in order 
to realise a faithful transfer of a many-qubit state, the relative phase of all 
qubits must be preserved. Therefore, the restrictive conditions of the perfect 
transfer do not suffice and it is necessary to hnd solutions that reconstruct 
both the amplitude and phase of a quantum state. To the best of our knowl¬ 
edge, general solutions have been derived only for 4- and 5-element chains 
(WGAs) with mirror symmetry [T7] . 

Here, we address the problem of the complete state reconstruction in hxed 
chains and hnd solutions by tailoring the coupling coefficients. We take the 
general approach to the inverse eigenvalue problem and analytically solve 
it for 4- and 5-element chains without symmetry constrains and for 7- and 
9-element chains symmetric around the centre. Unlike earlier approaches 
to the chain design that propose modihcation of the uniform chain [18] or 
solutions that asymptotically approach maximum transfer hdelity [19], the 
solutions proposed are exact and enable new design completely independent 
of uniform chains. We use the analytic solutions to numerically obtain the 
full state reconstruction in realistic WGAs realizable by direct laser writing 
[201 El]. Examples include the perfect transfer of a many-qubit state, an 
asymmetric array that supports the full state revivals and edge state Rabi 
coupling. Since the discretization inherently present in the analytical model 
approximates a WGA by a set of waveguides of inhnitesimally small diameter, 
we discuss limitations encountered in design of realistic WGAs. 

The paper is structured as follows. The mathematical model of a chain 
with next-neighbour linear coupling and condition of its commensurability 
are given in Section [^ In the same section, analytic solutions of the in¬ 
verse eigenvalue problem for general and symmetric chains are derived. In 
Section the analytic solutions are used to design realistic symmetric and 
asymmteric WGAs that support coherent transport with full state revivals. 
Further, an application of commensurable WGAs in construction of parallel 
data buses is suggested. Gonclusions are given and other chain conhgurations 
of interest are outlined in Section HI 
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2. Commensurable next-neighbour coupled chains 


Evolution of a state through a linear chain of n elements with the next- 
neighbour coupling is modelled by the Schrodinger equation of the form: 



( 1 ) 


where 'ip is the vector state in M” and An is a real nxn tridiagonal symmet¬ 
ric matrix that accounts for coupling between neighbouring elements. We 
assume that there are no loss or gain along the chain, hence that A is Her- 
mitian. In what follows, we will consider various n-dimensional cases of this 


eq.Q.” 


Periodicity of solutions requires all eigenvalues j = l,n of A to be 
commensurate, i.e. that there are integers rij such that 

T 

Vj : = (2) 

Period of oscillation is dehned as the smallest T for which the above equation 


holds for all eigenvalues of the chain. It is the smallest period with which 


the state can oscillate. However, if only some of the eigenvalues are excited, 
state revivals can occur with a longer period. 

Two and three element chains have unconditional periodicity. Their so¬ 
lutions are well-known and heavily applied in atomic and laser physics, see 
e.g. the theory of a three-level atom coupled by resonant laser light [22] • 

2.1. ^-element chain 

We use the shortest commensurate chain that permits incommensurate 
eigenvalues to illustrate the procedure used for solving the inverse eigenvalue 
problem. We start from a coupling matrix with real nearest-neighbour inter¬ 
action acting on of the general form: 


’ 0 O; 0 0 


^4(01, 02, 03) = 

0 og 0 oj 


O; 0 Og 0 


( 3 ) 


0 0 05 0 
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It is easy to verify that the eigenvalues a;i^2,3,4 of A4 are given by: 


^\/2a3^ + 2ai'^ + 2^/05'^ + 203^02^ — 203^0^^ + 

— ^\/2a52 + 20^2 + 2aj2 + 2^/a3^ + 2a3‘^a2‘^ — 2a3‘^ai‘^ + ag^ + _|_ 

|\/2a5^ + 2a2^ + 20^^ — 2\/ + 205^02^ — 2 a 3 ‘^ai‘^ + ag^ + 2 a 2 ‘^ai_‘^ + ai^ 

— \\/2a3'^ + 20^2 -|_ 2a_,2 — 2\fa3‘^ + 2a3’^a2^ — 2a3’^ai^ + + 2a2^ai^ + 

( 4 ) 

From ([^ we see that the chain is commensurate if and only if there are 
2 non-negative integers ni and n 2 with the greatest common divisor (GCD) 
equal to 1 for which 


\/2x^ A 2\/— Au^ = 27rni 

,_ ^ ( 5 ) 

\\/2x^ — 2\/x‘^ — Au'^ = 27rn2. 

In this case the eigenvalues are (±ni,±n2) scaled by an arbitrary scalar. 
Writing the frequency as g = ^ and solving the above equation we obtain 
that the coupling coefficients are 

ai = I (^ei \/ijii + n2Y - s"^ + e2\/ijii - n2Y - j 

02 = qs (6) 

03 = I (ei\/(w + ri2y - s2 - e2y/(ni - 77 , 2)2 - 

where ei and 62 in {— 1 , -(- 1 }, s an arbitrary real, and the frequency q = 2 tt/T 
an arbitrary non-zero real. 

Note that for any eigenspectrum given by ni and n2 an infinite number 
of chains can be constructed by choosing different s and 61,2- Not all values 
of s will give real values for the a^. However, all will give eigenvalues with 
the ratios {±ni,±n2}. 

An interesting example occurs when we choose degenerate eigenvalues, 
ni = n2- Then we are forced to choose s = 0 and hence 02 = 0 , meaning 
that the chain decomposes into two 2-element chains with the same coupling 
coefficients oi = 03. 
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2 . 2 . 5 -element ehain 

We turn to the commensurability of ^5(01, 02, 03, 04) given by 

’ 0 O; 0 0 O' 

Oi 0 Og 0 0 

^5(01,02,03,04)= 0 Og 0 05 0 . ( 7 ) 

0 0 % 0 

_ 0 0 0 0 _ 

The central eigenvalue of this chain in always 0 . We assume that the other 
eigenvalues relate to one another as (±ni, ±112) and repeat the procedure 
from the previous section. If af — 04 7^ 0 , then we obtain the commensurate 
array with the coupling coefficients in the form: 


„2 _ „2 

<^2 - Q y t‘2-s‘2 ) 

„2 _ 2 (-{nlnl-sH‘^)+{n\+nl-{s‘^+t^))t^\ 

03 “ y t2-s2 j 

O4 = q^t^ 

where s and t are arbitrary reals, and q is an arbitrary non-zero real. 
If — a| = 0 the solution assumed the form: 


al = q^n\ 



( 8 ) 


( 9 ) 


where q is an arbitrary non-zero real and is a rotation by an arbitrary 
angle (j). 

In the case of a chain that is mirror symmetric around the centre (a^ = 
ttn-j, j = 1, 2), the eigenvalues are easily calculated to be: {0, ±ai, ±A/2af+~af}. 
If Oi ■ 02 7^ 0 , the commensurability condition reduces to: 


aj = q^n\ 



( 10 ) 


If Oi = 0 , the system becomes a trivially periodic 3 -element chain with two 
outer elements remaining uncoupled to the others. If 02 = 0 the system 
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decomposes into two pairs of coupled elements that oscillate independently 
at both ends od the initially considered 5-element chain. 

For a chain with an equidistant energy spectrum, reverse engineering renders 
a solution from the family of Clebsch-Gordan coupling coefficients, e.g. for 

a symmetric 5-element array n 2 = 2ni, oi = gni, and 02 = Such 

chains describe coupling of atomic spins [8] and have been used to construct 
optical couplers im. 


2.3. Symmetric 7-element chain 

Now we turn our attention to the real, symmetric, nearest neighbour 
interaction with left-right symmetry acting on 


^47(01, 02 , 03) = 


0 

Oi 

0 

0 

0 

0 

0 


Oi 

0 

02 

0 

0 

0 

0 


0 

02 

0 

03 

0 

0 

0 


0 

0 

03 

0 

03 

0 

0 


0 

0 

0 

03 


0 

0 

0 

0 


0 02 
02 0 

0 Oi 


0 

0 

0 

0 

0 

QjI 

0 


( 11 ) 


and calculate the coupling coefficients that render the spectrum (0, ±ni, ± 77 , 2 , ± 77 - 3 ): 

If 03 7^ 0 : 

^2 


Of = q 




al = q^ 


nf — 


n^-hn^—nf 


( 12 ) 


o| = ,2 I 


where q is an arbitrary positive real. 

If 03 = 0, the chain is decomposed into two 3-element chains with period 

_ 27 r 

\/“ l +“2 

2 . 4 . Symmetric 9-element chain 


In the spirit of the previous sections, we solve the inverse problem for a 
mirror symmetric 9-element chain and obtain that the system with ^ 9 ( 01 , 02 , 03 , 04 ) 
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is commensurate if the Oj satisfy: 
If 03 7^ 0 and 04 7^ 0 : 


af = 




(n2+n^)(n§+n|-nf-n|)-(n|n|-nfn^) 


a| = 


n^n^— 71^712 


n^Ti2(n3+n|—n^—n^) 


( 7 i^^ 7 i 2 ){n^-\- 7 i^— 71 ^—7I2) —{71^71^—71^712) 


_ 2 f {‘nl+nl){nl+nl-nl-n^)-{nlnl-njn'^) 

^3-Q ^ni+ni-ni-ni) 


al = I n'^+nl-nl-n^ 


( 13 ) 


If 04 = 0 , the matrix consists of three diagonal blocks two of which are 4 - 
element chains and one block is the number 0. 

If 03 = 0 , there are also three diagonal blocks, that represent 3 -element chains. 


3. Coherent state transfer through optical waveguide arrays 

We apply the above theory to achieve coherent transport through optical 
WGAs. These arrays are accessible by modern direct laser writing [ 21 ] and 
lithography [ 23 ] fabrication techniques. Indeed, perfect transfer through the 
WGAs with equidistant eigenvalues has been demonstrated experimentally 
[mil]. The absence of loss due to the straightness of waveguides is one of 
important advantages of the WGA-based couplers over the circuits composed 
of nested optical couplers [ 21 ]. Here, we hrst describe the physical system 
in more detail and then apply exact general solutions to construct WGAs 
capable of full state revivals, a parallel data bus for optical quantum or 
classical computer, the hrst non-trivially commensurate asymmetric WGA 
and edge Rabi coupling. 

A typical dielectric optical waveguide consists of a region with refractive 
index higher than the index of the substrate that enables light conhnement in 
transversal direction and propagation along the guide. Waveguides couple to 
one another via evanescent helds with the strength that decays exponentially 
with the distance between them, [ 25 ]. Once the coupled coefficients of a 
chain are determined by the method described in Section]^ a WGA can be 
constructed by engineering the inter-waveguide separations as described in 

[min]- 
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Since the coupling coefficients are derived for a discrete model in Sec¬ 
tion 1^ that assumes a perfect transversal conhnement of a state to a chain 
element, we paid spatial attention to work close to this condition. This 
can be achieved in two ways: by working with high index-contrast waveg¬ 
uides and distant enough waveguides. Since the effective coupling in the 
former case requires extremely long propagation lengths, we chose to work 
with WGAs with large spacing between the guides. Fundamental modes of 
distant waveguides are strongly conhned to individual waveguides and the 
fundamental modes of closely-spaced waveguides form a supermode extend¬ 
ing through all of them. Fig. [^. All simulations with exception of the last 
example, were performed for the separation-to-diameter ratio around 2.2, 
which also secures weak waveguide coupling. 

The Rabi frequency q that features in expressions for matrix elements a* as 
a free scaling parameter is here set to 1. All simulated WGA were composed 
of identical circular waveguides with diameter 8.2 pm, core refractive index 
1.45 and substrate refractive index 1.445. Input modes were eigenmodes of 
individual waveguides at 1550 nm. To transfer the information through an 
array, we encode it in the amplitude and phase of a light wavefunction 'ip{z) = 

The time in Eq. [^is replaced by the 
space coordinate z along which the light propagates. Numerical simulations 
were performed using the hnite-difference beam propagation method and 
transparent boundary conditions f2U\ . 


3.1. Example 1: Phase and amplitude revivals 

An example of a symmetric 5-waveguide array inversely engineered to 
support full periodic state revivals is shown in hgure To prove the restora¬ 
tion of both the amplitude and phase of the wavefunction, we launch two 
input vector states of the same amplitude and different phases 0 and 7r/2 
into the WGA. The WGA is constructed to have eigenvalues 0, ±2, ±3 and 
to satisfy (10). The relative phase of input modes determines the propa¬ 
gation dynamics (light trajectory through a WGA). Nevertheless, the full 
revivals of the input state occur at dvrm, m = 1,2,3... irrespectively of the 
input phase, hgure The state vector redraws the same closed contour in 
the complex plane. The ’cleanness’ of the contour clearly distinguishes ar¬ 
rays with periodic dynamics from the arrays with quasi-periodic dynamics, 
whose state vectors never repeat the same path but hll in a subspace of the 
complex plane instead (see e.g. BZl). 
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Figure 1: (a) A waveguide array. Profiles of electric field amplitudes \E\ of supported 
modes sketched along z axis display evanescent field overlap responsible for coupling. 
Coupling coefficients ai correspond to those in matrices A. (b)-(d) Intensity (/ oc \E\^) 
profiles of the fundamental modes of pairs of circular waveguides separated by 2.7 D, 2.0 
D and 1.1 D, respectively, where D is the waveguide diameter. Red colour corresponds to 
the maximum and blue colour to the minimum intensity. White contours show waveguide 
cross sections. 


Further, use the transfer hdelity to estimate the accuracy required for re¬ 
vivals. It is dehned as F{z) = |('0(2;)|'^(O))— and is calculated for a range of 
input conditions. In the worst case,hdelity drops for 5% when the length of 
the array is 5% of the revival length. Taking into account that the usual dis¬ 
tance between the neighbouring waveguides is a few microns, a hdelity within 
1 % from 1 can be achieved with a resolution of a few tens of nanometers, 
which is achievable by the direct laser writing technique. 

3.2. Example 2: Parallel data bus 

The most exploited periodic solutions have been those used to realise 
the perfect state transfer from the beginning to the end of the chain. It 
implies the transfer with hdelity 1 of the state amplitude (energy) and is 
achieved if the ratio of diherences of subsequent eigenvalues ojn is rational, 
oji — Uj = (2m„ -b l)7r/Lo, where is an integer and Lq the chosen state 
transfer length m [2S]. The chain supporting such a transfer is referred 
to as quantum wire or quantum data bus. The procedure for the general 
inverse solution of this problem has been given in [13]. However, explicit 
solutions are still a few. Several systems have been proposed by choosing 
the eigenvalue scalings Un cc n, n^, n{n + 1) where n = 1, 2, 3.... The linear 
scaling renders a well-known array with spin-coupling coefficients formally 
equivalent to Clebsch-Gordan coefficients [29l[3n]. Analytical solutions have 
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wg1 wg2 wg3 wg4 wg5 wg1 wg2 wg3 wg4 wg5 




Figure 2: Results of numerical simulations of light propagation through a commensurable 
5-waveguide array with ni = 2, 712 = 3. Input states have different relative phases: (a,c) 
^(0) = (1,1, 0,0,0) and (b,d) tp{0) = (1, 0,0,0). (a,b) Upper graphs: intensity 

profiles (red colour corresponds to the highest intensity). Lower graph: |'i/’(-z)P obtained 
by solving Q. (c,d) Rotation of the state vector ip in complex plane. Colour coding 
of the lines is blue, green, red, cyan and purple going from the bottom waveguide up, 
respectively. 


been reported for short chains with 6, 7, 8 and 16 elements [3T]. In laboratory, 
the perfect state transfer through chains with more than 3 elements has been 
demonstrated in optical waveguide arrays (WGAs) [Ij, and between Zeeman 
states of atoms constituting a Bose-Einstein condensate [S]. Realisations in 
other physical systems, such as quantum dots, have been proposed [32] . 

Here we use the commensurability conditions to expand the set of known 
analytic solutions for perfect transfer. Since a WGA can support periodic 
dynamics during which the perfect transfer to the end waveguide never oc¬ 
curs, the condition of perfect transfer is not equivalent to the periodicity 


11 









































condition and the explicit solutions for perfect transfer given below are sub¬ 
sets of solutions derived in Section The perfect transfer can be realized in 
S-elements WGAs that satisfy ( 10 ) and have even 77-2 and ni < 77,2. Similarly, 
the condition ( 12 ) and even tti, odd 772, odd 773, and 772 < 77 i < 773, guarantee 
the perfect transfer between the end waveguides of a 7 -element array. These 
solutions include known Heisenberg chains with equidistant eigenvalues. An 
example of such a chain with (771 = 2 , 772 = 1 , 773 = 3 ) is shown in figure]^. 

Further, we propose a parallel optical data bus that can transfer a number 
of states simultaneously. It is based on a commensurate WGA acting as a 
multiport device that supports both the parallel input and read-out from 
77 waveguides. Here, instead of limiting the consideration to the perfect 
(energy) transfer from one to the other end of a WGA, we exploit coherent 
transfer of the complex wavefunction 'ip for all possible input states. Any 
WGA with a commensurate eigenspectrum can perform this function. In 
Fig-i we show two examples. The hrst is the Heisenberg-like chain and 
the second is an array with Fibonacci sequence of eigenvalues 773 = 13 , 772 = 
5 , 773 = 8. Of course, the eigenvalues do not have to posses any particular 
regularity apart from that dictated by the formulas in Section and the 
condition of no loss/gain. Indeed, all WGAs given as examples throughout 
this Section are potential solutions for construction of a parallel bus. 




Figure 3: Results of numerical simulations showing the complete transfer through a 7- 
waveguide array with (a) equidistant eigenvalues ni = 1, 772 = 2, 773 = 3 and input vector 
state ■!/>(0) = (1,0,0,0,1,0,0), (b) Fibonacci eigenvalues rii = 3, 712 = 8, 773 = 5. The 
input vector state is (1,0, 0, 0, 0,0,0). \'ip{z)\^ obtained by solving Q. Colour coding of 
the lines is blue, green, red, cyan, purple, olive green and black going from the bottom 
waveguide up, respectively. 
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Another chain that supports a parallel read-out are atoms with different 
spins that can be imaged simultaneously using the Stern-Gerlach scheme [8]. 
However, construction of an arbitrary input state comprising different atomic 
spin levels requires sophisticated coherent control techniques, [7]. 

3.3. Example 3: Asymmetric arrays 

While the chains reported in literature treat only symmetric arrays, we 
observe that the full state revivals are possible in asymmetric arrays that 
are constructed without any a priori restrictions on coupling coefficients. 
Solutions in Section m and |2.2| enable construction of such arrays. An 
example in Fig. shows a 4-guide array with commensurate eigenvalues 
Hi = ±1, 77.3 = ±3, and arbitrarily chosen parameters s = \/3, ^1 = 1 
and €2 = —1. It is composed of the WGs separated by L 2,3 = 0.977Li 2 and 
-^ 3,4 = O. 954 L 1 2 . To produce the required small difference in inter-waveguide 
separations, a resolution of a few percent is required. For usually used sepa¬ 
rations of several microns, such resolution is achievable by femtosecond laser 
fabrication technique. 



Figure 4: Results of numerical simulations of light propagation through a commensurable 
asymmetric 4-waveguide array with eigenvalues set by ni = 1, n 2 = 3 and the input vector 
state tp{0) = (1, 0,0,0). Upper graph shows intensity prohles (red colour corresponds to the 
highest intensity). Waveguides are shown by black lines. Lower graph: obtained 

by solving Q. Colour coding of the lines is blue, green, red and cyan from the bottom 
waveguide up, respectively. Small deviations from the periodic propagation are due to the 
limited mode confinement in a realistic waveguide and numerical artifacts. 
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3.4- Example 4- Reduced arrays 


We further explore asymptotic cases of effective reduction of an array in 
terms of the number of elements. When some inter-waveguide separations are 
much smaller than the others, propagation dynamics comprises slow and fast 
oscillating helds, favourizing propagation through the closely spaced waveg¬ 
uides. For instance, in a 7-waveguide array with eigenvalues that satisfy 
the condition » rii » 77 - 2 , the light power is conhned to the 4 outer 
waveguides at any point of propagation, hgure According to the formulas 
given in the above condition renders an array with a stronger coupling 


between the inner than the outer elements, 03 >> 02 >> Oi, with the den¬ 
sity of waveguides decreasing from the centre of the array outwards. Note 
that, if cut at the any length that corresponds to an odd multiple of 7r/2, 
this WGA can be used as an equal 1x4 splitter (splits a signal into 4 with 
equal amplitudes). 


4. Conclusions 

We have investigated hnite chains with linear coupling as means of coher¬ 
ent state transfer in quantum computers. Coherent propagation with peri¬ 
odic state revivals is achieved by reversely engineering couplings between the 
neighbouring elements of a chain to yield its eigenspectrum commensurate. 
The inverse eigenvalue problem is solved analytically for general chains with 
4 and 5 elements and for symmetric chains with 7 and 9 elements. These so¬ 
lutions comprise particular cases actively studied to achieve the perfect state 
transfer from one to the other end of the chain. However, the presented solu¬ 
tions provide a much wider and still under-explored possibilities for coherent 
propagation and revivals of both the state amplitude and phase through a 
set of parallel channels. The full inverse solution is demonstrated on experi¬ 
mentally accessible optical WGAs. In such a system, coupling coefficients are 
controlled by tailoring the inter-waveguide distances. Practical limitations 
neglected in a discrete analytical model are considered. Presented analyti¬ 
cal and numerical solutions offer a plenitude of possibilities for engineering of 
WGAs that are now realizable by mature laser fabrication technologies. They 
are characterised by low loss and robustness of architecture. Finally, besides 
the mathematical challenge to hnd analytical solutions for linear chains with 
larger and even number of elements, there is a considerable interest in co¬ 
herent transfer through two dimensional and closed chains, e.g. multicore 
fibres. 
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Figure 5: Results of numerical simulations showing light dynamics in 7-waveguide array 
with ni = 100, 712 = 1, = 10000 that reduces to a 4-waveguide array. The input vector 

state is (1,0,0,0, 0, 0,0). obtained by solving Q. Colour coding of the lines is 

blue, green, red, cyan, purple, olive green and black going from the bottom waveguide up, 
respectively. 
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